Electrostatic screening in fullerene molecules. 
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Abstract 



We study the Coulomb interactions in fullerene molecules within a con- 
tinuum formalism. The model gives rise to a renormalizable field theory, 
which has many similarities to standard quantum electrodynamics. The ef- 
fective electric charge at low energies is reduced by screening processes. The 
associated renormalization of the one electron Green's function leads to the 
vanishing of the quasiparticle pole. It implies the dissappearance of coherent 
one particle excitations, in close analogy to the one dimensional Luttinger liq- 
uid. The relevance of these results for Ceo and related molecules is discussed. 

75.10.Jm, 75.10.Lp, 75.30.Ds. 
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I. INTRODUCTION. 



The localization of the electrons within the Cgo spheres implies that electrostatic inter- 
actions cannot be neglected in fullerene crystals. The simplest estimate for the charging 
energy of a C 60 molecule, e 2 /R, where R is the radius, gives a rather large value, ~ 4eV. 
Moreover, the low density of states at the Fermi level in graphite reduces metallic screening 
(see below). In general, the intermediate size of the molecule puts it in a regime different 
from that for an isolated atom or a bulk system. Different approaches have been used in the 
study of electrostatic effects It has also been proposed that they can be the origin of 

superconductivity in doped systems 

In the present work, we use a simple model for the long wavelength electronic properties 
of fullerene molecules [PH|,PJ], to analyze the screening properties. It has been shown that it 



gives a reasonable approximation to the electronic levels of Cgo- Furthermore, the scheme 
is sufficiently general to describe other systems with similar topology, like nanotubes and 
fullerene 'onions'. The model can also be used to study other features, like the electron- 



phonon interaction with reasonable accuracy, a shown in a previous work [12]. 

Calculations based on this model for the electrostatic interactions between the highest 
occupied orbitals (in doped systems) are in good agreement with more detailed calculations 
i- ' 

In the next section, we discuss how to incorporate the Coulomb interaction into this 
model. It will be shown to lead to a renormalizable field theory in (2+1) dimensions. Then, 
we analyze in detail the leading divergent diagrams in the theory, and show how to deal 
with these divergences by definining appropiate renormalized couplings. The flow of these 
couplings, and its implication for some quantities of physical interest, is discussed in the next 
section. Finally, in the conclusions, we discuss the applicability of our results to C 60 and 
related systems, as well as the changes induced by the finite curvature and doping expected 
in these systems. 



II. THE MODEL. 

The low lying electronic states of an isolated graphitic sheet are well approximated by 
an effective Dirac equation in (2+1) dimensions, using Huckel's theory of conjugated carbon 
compounds. This description plays the same role as the effective mass theories used for the 
study of the long wavelength properties in other materials. The existence of rings with an odd 
number of atoms in Cqq leads to additional complications, which require the introduction of 
a fictitious gauge field, to account for the frustration that these rings induce in the electronic 
states. 

We use the preceding scheme to analyze first the screening properties of an isolated 
graphite plane. The single sheet, in the absence of interactions, is semimetallic, with a 
vanishing density of states at the Fermi level. Hence, the screening properties are expected 
to differ significantly from those typical of standard metals, as no screening length, kpT — 
y / 47re 2 A^(ei?), can be defined. 

Within the approximation described above, we can write the total hamiltonian, including 
the Coulomb interaction, as: 
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2 J J 4tt | ri — r 2 | V ; 

where the first part describes the independent electron hamiltonian by means of the Dirac 
equation. The 7's are Dirac matrices, t is the hopping between nearest neighbor orbitals, 
~ 2.2eV and a is the interatomic distance, ~ 1.4A. The Fermi velocity is given by 3ta/(2h). 
The second part is the standard electrostatic interaction. 

We can now try to analyze the interaction term perturbatively. To lowest order, we have 
the exchange term, depicted in figure (1). Note that what we define as Coulomb interaction 
is not the same as what is traditionally assumed for 2D charges, which is a logarithmic 
potential. We suppose that the electronic charge, although localized within a plane, gives 
rise to a 3D electrostatic potential, which decays as r _1 . This difference is crucial to the cutoff 
dependence of the diagrams that we will encounter. For instance, the exchange contribution 
to the self energy, figure (1), can be writen as: 

-> -» r ulc cPk 1 

E xc (k) oc e 2 7 • kj —3- (2) 

where u c is an upper cutoff, above which the linearization we use ceases to be valid. The 
integral in (2) diverges logarithmically as function of its upper and lower limits. 

This divergence is an indication of the need of some kind of renormalization procedure. 
The possibility of logarithmic corrections to the physical quantities can also be inferred from 
standard power counting. If we choose units such that velocities are dimensionless, then the 
coupling constant, e 2 , is also dimensionless. In fact, this power counting is identical to that 
in standard quantum electrodynamics. 

In the following, we will account for the divergence in equation (2), and those other 
divergences which appear at higher orders in e 2 , by renormalizing the terms in the initial 
hamiltonian. A proof that it can be done, and a detailed explanation of most technical 
steps, can be found in Thus, many intermediate steps will be omitted in this paper, 



athough we will discuss in full the physical implications. 

In order to simplify the calculations of the divergences, we will couple the 2D Dirac 
equation in our theory to 3D photons, which then induce the Coulomb interaction, eq. (1). 
This step is not simply a technicality, as it also allows us to discuss retardation effects. 
Retardation due to the finiteness of the velocity of light is relevant at sufficiently large 
distances. This is a well known effect in the study of Van der Waals forces, for instance. If 
the typical separations, I, are of the order of Ae/ (he), where Ae is a typical energy splitting, 
and c the velocity of light, then the dependence of the interaction on I changes. We study 
the generalization of these interactions to graphite planes, and, as we will see, similar effects 
appear. 

To include the 3D electromagnetic field, we generalize the hamiltonian, eq. (1), to: 

H=haJ d 2 r ^(r) 7 • Vtf(r) - e J d 2 r (3) 

where we are using units such the velocity of light, c = 1, and we are omitting the 
term associated to the free electromagnetic field. It is important to realize that this model, 
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although formally similar to QED, does not exhibit Lorentz invariance. We have two dimen- 
sionless velocities, Vp = 3ta/(2h), and c. Initially, vp/c <C 1, which corresponds to the non 
relativistic limit. 

The interaction of the field and the electrons in the layer is described in the standard 
fashion, by coupling to the conserved current: 

j„ ~ (i ^70^, i u*7*) (4) 

This poses some technical problems since the electromagnetic field propagates in three- 
dimensional space while we want the dynamics of the electrons to be confined to the two- 
dimensional layer. Although this may not be achieved in general, it turns out to be possible 
by specializing to the Feynman gauge, which enforces the constraint: 

= (5) 

In order to proceed perturbatively, we need to know the propagator of the field in 
this gauge, which is: 

, pik-(r-r') -iu(t-t') 
(TA At , r) M t>, r<)> = -«„ / ——^ r -— (6) 

which then couples to the 2 + 1 dimensional current. The nonrelativistic limit of any 
quantity is obtained by expanding in powers of Vp/c. As mentioned above, we may ex- 
pect differences between treatments based in hamiltonian (1) and hamiltonian (3), when 
retardation effects become relevant. 



III. CALCULATION OF THE LEADING DIVERGENCES. 

We expect logarithmic divergences to all orders in a perturbation expansion in powers 
of e 2 . To facilitate the computation, we use a dimensional regularization procedure [|13|,[n[ . 
The bare parameters in the hamiltonian undergo multiplicative renormalization: 

Vbare = Z v Vr 
&bare Z e &R 
Vbare = Z^ R 

Kare = Zf A» R (7) 

A detailed calculation of the leading divergences in each of these factors can be found in 
13fl . As an example, to lowest order in e 2 , we find the following contributions: 



^ = l + ^(l-2^ 2 ) f l dx . VT ~^ 2 - + 0(e 
87r^ v ' Jo 1 — x + v z x e 

Z„ = 1 / dx- 



7T 2 Jo (1 — x + v 2 x) 2 e 
; 2 u n _ 2X f 1 , VT^x 1 



l-2v 2 ) / dx-^ — - + 0(e 4 ) (8) 



87r 2 \ 'Jo 1 — x + v 2 x e 
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where v stands for the Fermi velocity in dimensions where c—1. In the nonrelativistic 
limit, v — > 0, we can write a an expansion in terms of v. The leading term in Z v is: 

Z v = l- -^_e 2 L l -F ( 1 I v 2 ) - 2nv (l - 2v") F (-, -; 2; W 2 ) 
16tt 2 I v V2'2'2' / V y V2'2'2' / 

+4 (l - 2t; 2 ) F (l, 1; 1; /) - 4F (l, 2; |^ 2 ) } 1 + O (e 4 ) (9) 
This expression is written in terms of hyp ergeo metric functions. 

The model exhibits gauge symmetry, which is not broken by our choice of regularization 
scheme. That imposes constraints on the values of the Z's. The renormalized Lagrangian 
can be written as: 

Share = Z kin J dtSr ^(-7 <9 + Z v V^f • V)* 

-Z int ie J dtd 2 r ^(-70^0 + Z v v*y • A)V (10) 

— 1/2 

and the gauge invariance implies that Zkin = Zi nt . This, in turn, implies that Z e = Z A . 
Furthermore, we can infer from the structure of vortex diagrams, figure (2), that the electric 
charge is not renormalized. The reason is the insertion of a 7 oerator at each vortex, which 
makes the diagram tend to zero as the momentum transferred, q — * 0. We have checked 
that these requirements are satisfied to order e 4 . 



IV. RENORMALIZATION FLOW. 

Once the multiplicative constants, eq. (8) are known, we can analyze the dependence of 
the propagators upon a change in the cutoff: 

U A + e 2 )|- + p e (v, e 2 )^ - 7 („, e 2 )) G(u, k, A; v, e 2 ) = (11) 



where 



(3 v (v,e 2 ) = A ^ V R 

1 (v,e 2 ) = A-^logZ^ (12) 



The well-known solution to this equation takes the form 



G(u, k, pA; v, e 2 ) = exp y l^j G(u, k, A; v eff (p), e 2 eff (p)) (13) 



and the effective parameters are given by 
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d 

p-Q- p V eff(p) = -Pv(v eff ,4 ff ) 

d 

Pfy e eff{p) = -Pe(v ef f, e 2 eff ) (14) 



Increasing the scale p of the cutoff is equivalent to measure the observables of the theory 
at large distance scale. From the one-loop order results of the previous section we see that 
e e ff remains constant, while v e ff does not. We expect, therefore, the effective coupling of 
the electronic interaction e 2 /(4nv) to have a nontrivial renormalization group flow in the 
infrared regime. Actually, from equations 7 and 12, 
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8 5^ 1 , _j_ , 

dp v R lQnv eff 




so that the asymptotic behaviour of the coupling is: 

1 e 2 1 e 2 / 1 e 2 , A x ~ 1/2 



(A/A ) = 1 log— (16) 

Here, A is the energy at which the effective parameters are defined. 

The behavior of the coupling constants as function of energy, eq. fll6|), is a consequence 
of the screening processes which take place within the graphite sheet. The effective coupling 
is reduced at low energies and long distances. The increase in the Fermi velocity leads to 
a reduction in the density of states near the Fermi points, which can be thought of as the 
tendency towards the formation of a Coulomb gap. A true gap, however, does not exist. 

The fact that the electromagnetic interaction becomes comparatively weaker at long 
distances implies that our results, although perturbative, are qualitatively correct at low 
energies. The initial coupling is large in dimensionless units (see below), but is reduced 
as the renormalization proceeds. Hence, below a given scale, our results, which have been 
calculated to order e 4 , are valid. We can not predict quantitavely, however, the relation 
between the coupling at that scale and its bare value. 

The flow of the parameters towards the weak coupling fixed point resembles the flow 
in QED. In fact, even a "Landau pole" appears, if the equations are iterated towards the 
ultraviolet. This pole is not too far above the bare parameters. It probably signals a phase 
transition for sufficiently large initial couplings. We should also remark that the flow of the 
Fermi velocity towards large values is bounded by the velocity of light. At this point, the 
model becomes Lorentz invariant. This symmetry should be preserved by the subsequent 
renormalization flow. 

The fact that the propagators acquire non trivial exponents should be reflected in all 
physical properties. For instance, the temperature dependence of the specific heat will differ 
from the T 2 law expected in this semimetal, and the exponent will no longer be an integer. 

The one electron propagator is of special importance, as it measures the coherence of 
the quasiparticles, and also governs the hopping processes to other molecules or planes away 
from the one being considered. The local density of states in the non interacting case is 
lmG(u) ~ |co>|. This power law is changed to lmG(uj) ~ |cj| 1+ °( e '( ta ) \ This appearance of 
anomalous exponents resembles closely the features of a one dimensional Luttinger liquid. 
It reflects the same physical processes, namely, the absence of a coherent quasiparticle pole 
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near the Fermi level. Thus, the system is not truly metallic, although it does not exhibit a 
gap either. This result implies, among other things, a significant reduction in the effective 
hopping elements to an external system, as discussed extensively in connection to Luttinger 
liquids [|L6[| . 



V. APPLICATION TO FULLERENE MOLECULES. 

In fullerene molecules, a number of factors not considered so far play a role: 

- The initial value of the coupling constant, e 2 /(ta) ~ 4 is large, although in a per- 
turbative scheme it appears divided by numerical factors greater than unity. As mentioned 
earlier, the existence of a large initial coupling does not change qualitatively the low energy 
features of the model. Screening processes reduce this values, until a scale is reached be- 
low which perturbation theory is applicable. The actual exponents associated to physical 
quantities, however, will depend on the value of the coupling at this scale, which cannot be 
obtained perturbatively. 

- In a finite system, the typical dimensions act like an infrared cutoff, below which 
scaling is not applicable. This is a severe restriction for the Ceo molecule, whose radius, 
R ~ 3.2Ais close to the distance between carbon atoms, which is the initial cutoff. Thus, the 
renormalization scheme is only a qualitative guide to the role of electron-electron interactions 
in this molecule. It suffices, however, to rule out models based on the existence of strong 
metallic screening beyond the C-C separation, like the Hubbard model. The renormalization 
flow, on the other hand, should be more precise for the larger fullerene molecules, which 
may exhibit the unconventional features reported here. 

- In doped systems, the finite density of states at the Fermi level allows for metallic 
screening. The best known cases are the family of alkali doped Ceo, A3C60, where A stands 
for an alkali metal. Then, the doping concentration is 1/40 electrons per site and per spin. 



That means that the Fermi is equal to k F = a 1 \J 4^/(15-^/3) ~ 2A. Using the density of 
states for a linear dispersion relation, we can define a screening length: kp% ~ 3.4 A. This 
number is similar than the radius of the molecule, so that metallic screening should not 
play a significant role. Alternatively, we can use the density of states of K 3 C 60 from band 
structure calculations ||15|| . Its value is 25 eV _1 per Cqo molecule. This number is higher 
than the estimate given above. Hence, we find kpp ~ lA. However, as discussed in [|15||, the 



high density of states is due to the hybridization of a single orbital from each Ceo molecule. 
The charge transfer between Ceo molecules in the presence of an electric field does not imply 
the existence of local, intramolecular rearrangements of the charge distribution, because the 
charge remains in this triply degenerate orbital. 

These points restrict the application our results, obtained for infinite graphite sheets, to 
the Ceo compounds. However, some qualitative conclusions remain valid: 

- The long range effects of the Coulomb interaction are significant, and cannot be 
neglected. From the preceeding sections, we can conclude that there is no screening length 
in isolated sheets of graphite. The inclusion of finite size effects, and of the doping does not 
change this result. The new length scales which appear when these phenomena are taken 
into account are of the order of the Ceo radius. Thus, the effective range of the Coulomb 
interaction is, at least, comparable to the size of each Ceo unit. 
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- The absence of coherent quasiparticles implies the existence of strong shake up 
effects when electrons tunnel into the molecule. The behavior of the Cqq spheres resemble 
that of finite ID Luttinger liquids. As discussed above, the flow of the coupling constants is 
towards a fixed point with no quasiparticle pole. Coherent quasiparticles cannot be defined 
at this fixed point. Finite size effects halt this flow before this situation is reached. However, 
its effects will be manifested in a variety of physical processes. The most relevant of them is, 
possibly, the reduction of the effective tunneling matrix between neighboring spheres. This 
is a consequence of the strong shake-up effects which are induced when a single electron is 
transferred between spheres JT7|-P^[. The ration between the bare intermolecular tunneling 
element, t inter , and the effective tunneling at low energies, t e ff, goes as t e ff/t in t er ~ (u/uj c ) 9 , 
where u c is the initial ultraviolet cutoff, proportional to the it bandwidth, and g is the 
exponent which characterizes the anomalous scaling of the density of states at low energies. 
As mentioned earlier, this scaling low cannot be continued to energies below a given cutoff. 
This cutoff is of the same order of magnitude as ta/R, where t is the hopping between 
nearest neighbor tt orbitals on the same sphere, a the interatomic distance, and R the radius 
of the molecule. Thus, u;/u; c ~a/i?~ 1/3. From the fact that we are in the strong coupling 
regime, e 2 /(at) > 1, we infer that g is of order unity. Hence, these effects play a significant 
effect in renormalizing the effective hopping between Cqq molecules. Note that shake up 
effects cannot be studied within conventional band structure calculations. 
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FIGURES 

Figure 1. Lowest order correction to the electron propagator (exchange term) 
Figure 2. Lowest order correction to the electronic charge. 
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